Infrared propagators of Yang-Mills theory from perturbation theory 
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We show that the correlation functions of ghosts and gluons for the pure Yang-Mills theory in 
Landau gauge can be accurately reproduced for all momenta by a one-loop calculation. The key point 
is to use a massive extension of the Faddeev-Popov action. The agreement with lattice simulation 
is excellent in d = 4. The one-loop calculation also reproduces all the characteristic features of the 
lattice simulations in d — 3 and naturally explains the pecularities of the propagators in d — 2. 
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The infrared (IR) physics of strong interaction is well 
described today by lattice simulations of Quantum Chro- 
modynamics (QCD). This tool is now commonly used to 
determine the spectrum of particles, cross sections and 
other physical observables (see for example [l|). The an- 
alytical (or semi-analytical) approaches have not reached 
such a high level of development, mainly because the 
standard perturbative approach breaks down at low ener- 
gies. This calls for more sophisticated techniques, and a 
good guiding principle in developing them is to compare 
their predictions with lattice simulations. Unfortunately, 
the simplest quantities that can be computed in ana- 
lytical approaches are non-gauge-invariant and therefore 
require to fix the gauge. For this reason, a considerable 
amount of work has been performed in the last few years 
to study gauge-fixed versions of QCD in the lattice. 

The simplest quantities that can be analytically stud- 
ied are the 2-point correlation functions in Landau gauge, 
and most of the gauge-fixed simulations have focused on 
these quantities. In the case of pure gluodynamics (with 
no quarks), some facts are now clearly established. First, 
the gluon propagator does not diverge in the IR but tends 
to a positive constant for d > 2 [2j-|4| and to zero in 
d = 2 [a, E| . The ghost propagator is divergent in the IR 
limit with an enhancement when compared to the free 
propagator: the dressing function (i.e. the propagator 
times momentum squared) is monotonically decreasing 
with momentum. It seems to approach a finite posi- 
tive constant in the IR for d > 2 and diverges in this 
limit for d = 2. It is also well documented that the 
Kallen-Lhemann spectral function associated with the 
gluon propagator is not definite positive @, Q . 

Let us now recall the various analytical approaches 
that have been used to determine these correlation func- 
tions. The standard perturbation theory in the frame- 
work of Faddeev-Popov (FP) gauge-fixing, as is well- 
known, is unable to access the IR limit of the theory 
because it presents a Landau pole. This may be related 
to the fact that the FP procedure does not fix completely 
the gauge because of the Gribov ambiguity: there ex- 
ists in general several gauge transformed configurations 
(Gribov copies) that satisfy a gauge condition [9|. A line 



of investigation has been developed to restrict the func- 
tional integral in order to take into account only a subset 
of the Gribov copies (hopefully, only one). This leads 
to the Gribov-Zwanziger model [9l4llj and some variants 
of it [12j . The IR propagators have also been studied by 
using Schwinger-Dyson and Non-Perturbative Rcnormal- 
ization Group equations. In these approaches, one solves 
a truncated version of an infinite set of coupled equations 
for the vertex functions. Depending on how one imple- 
ments these ideas, two families of solutions have been 
found: i) the scaling solution ljj |l4|, where the gluon 
propagator goes to zero in all dimensions and the ghost 
propagator is more singular than the bare one in the IR, 
and ii) the decoupling solution [l5| where the propaga- 
tors have behaviors in qualitative agreements with the 
lattice simulations. We note at this level that all these 
approaches lead to quite involved calculations, with in 
some cases an important numerical part. 

In this letter we take a more pragmatic point of view. 
We do not try to find a gauge-fixed theory that would 
be justified from first principles, but propose a minimal 
modification of the FP action that can account for the 
lattice simulation results. Of course, this phenomeno- 
logical approach can only be motivated a-posteriori, if 
it describes in a satisfactory way the simulation results. 
Our main guide is the observation that the gluon propa- 
gator tends to a finite positive value in the IR for d > 2. 
We propose to impose this property at the tree level by 
adding a mass term for the gluon in the FP action }25| . 
We do not change the ghost sector since the ghost propa- 
gator is found to be IR divergent in the simulations. This 
leads us to consider the Landau-gauge FP euclidean La- 
grangian supplemented with a gluon mass term: 
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where {D^Y = <V° + g$ ahc A\c c and the field strength 
F£„ = d M A v - d„A„ + gf abc A b ^A c u are expressed in 
terms of the coupling constant g. The Lagrangian (JTJ 
corresponds to a particular case of the Curci-Ferrari 
model [l6|. At the tree level, the gluon propagator is 



2 



massive and transverse in momentum space: 

1 



P + Ttl 



(2) 



with Pi, (p) 



-p^p v /p 2 . It is interesting to note that 



the spectral density associated with the propagator ([2]) is 
positive and therefore there is no violation of positivity 
at the tree level. We conclude that violation of positivity, 
if it exists in this model, is caused by fluctuations. 

Actually, the gluon propagator observed in the lattice 
is not compatible with the bare propagator ([2]) and we 
will show below that, by including the one-loop correc- 
tions, one obtains propagators for gluons and ghosts that 
are in impressive agreement with those obtained in the 
lattice in d — 4 (including positivity violations) and that 
reproduce at the qualitative level the results for d = 3. 
Let us mention that a mass term has been used to im- 
prove perturbative QCD results in order to reproduce 
the phenomenology of Strong Interactions 17[. More- 
over, there are successful confinement models I8[ that 
use actions including a gluon mass term. 

When analyzing the model described above, we must 
face the problem that the mass term breaks the BRST 
symmetry [l9| which is very important in the perturba- 
tive analysis. This symmetry has the form 



SA% = r,(D li c) a , 5c a = - v Zf 



abc^b 



Sc a — T]ih a , 



6ih a = 0, 



(3) 



where 77 is a global grassmanian parameter. The BRST 
symmetry is in general used to prove the renormalizabil- 
ity of the theory. However, the breaking of the BRST 
symmetry by the mass term is soft and therefore does not 
spoil renormalizability [lij . The BRST symmetry is also 
used to reduce the state space to the physical space, in 
which the theory is unitary (at least at the perturbative 
level) and the breaking this symmetry spoils the stan- 
dard proof of unitarity. This problem is actually com- 
mon to essentially all methods that try to go beyond the 
standard perturbation theory (as the Gribov-Zwanziger 
model) because they all break the standard BRST sym- 
metry. In this respect, the model considered here is not 
in a worse position than other approaches considered in 
the field. We must stress that this model is equivalent to 
the standard FP model in the ultraviolet limit p » Aqc d 
if m ~ Aqcd . This means that in the domain of validity 
of standard perturbation theory, the model is as unitary 
as QCD. The unitarity of the model in other momentum 
regimes is of course an important open problem, as it is 
in all gauge fixings in which standard BRST symmetry 
is broken. 

The model with Lagrangian ([1]), as a particular case of 
the Curci-Ferrari model, has a pseudo-BRST symmetry 
(not nilpotent) that has the same form as the standard 
BRST © except for the h variation which reads 8ih a = 
rjm 2 c a . On top of this symmetry, the Lagrangian has all 
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FIG. 1: Four-dimensional correlation functions for SU(2) 
gauge group. The results of the present work (red curve) are 
compared with lattice data of (green points). Top figure: 
gluon propagator. Bottom figure: ghost dressing function. 
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FIG. 2: Four-dimensional gluon propagator for SU(3) gauge 
group times p 2 . The results of the present work (red curve) are 
compared with lattice data of Q (green open circles) and 0] 
(green crosses). The (black) dashed curve is the ultraviolet- 
improved curve obtained by the renormalization group. 



the standard symmetries of the FP action for the Landau 
gauge. This includes the shift in antighost c — > c + cst., a 
symplectic gro up l20i| -, and four gauged supersymmetries 
recently found [21| . As a consequence, the mass [22j and 
coupling constant 23[ renormalization factors (even in 
presence of the mass term 2l| ) are fixed in terms of gluon 
and ghost field renormalizations. 

We present now the 1-loop calculation of the propaga- 
tors, which requires the calculation of four Feynman dia- 



grams. It is convenient to parametrize the gluon G'^Kp) 
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and the ghost G (p) propagators in the form: 

G ab (p) = S ab F(p)/p 2 , G% (p) = ( P )5 ab G(p). (4) 

The F(p) is known as the ghost dressing function and the 
scalar function G{p) will be refered to as the gluon prop- 
agator below. We choose the following renormalization 
conditions: 



G(P = 0) : 

F(p = ,i) 



1/to 2 , 
1. 



G(p = /i) = l/(m 2 + /x 2 ), 



(5) 



We use the gluon-ghost vertex in the Taylor scheme [23| 
for the coupling constant g. 

We consider first the 4-dimensional case. The one- 
loop result for the renormalized functions F(p) and G(p) 
(imposing the renormalization prescriptions ([5])) are: 

G- 1 (p)/m 2 = S + l + |^{lll S - 1 -2 S - 2 

+ (2 - s 2 ) log s + 2(s" 1 + l) 3 (s 2 - 10s + 1) log(l + s) 

+ (4s- 1 + I) 3 / 2 (s 2 - 20s + 12) log f ^±I~^ ) 

V V 4 + s + V« / 

-(s^/to 2 )} 

^(p) = 1 + - slogs + (s + l) 3 s- 2 log(s + 1) 

- s- 1 - (s -> /; 2 /to 2 )} (6) 

where s = p 2 /m 2 . 

In Fig. [U we compare these expressions for the SU(2) 
gauge group with the lattice simulations of ;2] . The best 
choice of parameter is g — 7.5 and m — 0.68 GeV when 
normalization prescriptions are imposed at /J, = 1 GeV. 
One observes that both gluon and ghost propagators can 
be fitted with the same choice of parameters in a very 
satisfactory way. Note that the normalization conditions 
of the lattice simulations are not compatible with ([5]) so 
that we have to introduce a global multiplicative renor- 
malization factor when comparing the curves. 

We have also compared our results with the data of 
two different lattice studies d i] for the SU(3) group. 
The two data sets have different overall momentum scale 
and we have rescaled the momenta of the data of [H 
for superimposing them with those of We represent 
in Fig. [2] the dressing function of the gluon instead of 
the propagator in order to make visible the ultraviolet 
regime. The best choice of parameters is g = 4.9 and 
to = 0.54 GeV (again with /1 = 1 GeV) and it leads to 
a very satisfying agreement for momenta p < 2 GeV. It 
is important to stress that expressions (|6]) are 1-loop re- 
sults obtained from a fixed coupling constant calculation 
in a fixed renormalization point. It is well-known that 
in order to analyze the regime p 3> to, one must take 
into account renormalization group effects and in par- 
ticular the running of the coupling. The corresponding 
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FIG. 3: Four-dimensional real space propagator C(t) for 
SU(3) gauge group. The curves grows when t tends to zero, 
saturating at the value C(0) ~ 2.1. 



procedure is standard and once it is implemented (see 
[24| for details) , the agreement is essentially within error 
bars for p > m as is also shown in Fig. [2] A very good 
agreement is also obtained for the ghost dressing func- 
tion In any case, it is obvious that when p ^ to, 
the model ([1]) reproduces correctly the high momentum 
regime once renormalization group effects are taken into 
account. 

An interesting feature of the 1-loop gluon propaga- 
tor is that it is non- monotonous in the IR. In fact, 
the inverse propagator behaves at small momenta as 
to 2 + V 5 V/(1927r 2 )log(p 2 /TO 2 ) + 0(p 2 ). This predic- 
tion of our calculation is very small for d = 4 and it is 
not visible in Fig. [1] but appears clearly in d = 3, see 
below. 

An important property of the propagators measured 
on the lattice is the violation of positivity. One way to 
extract it is to calculate the quantity: 



C(t) 



dp 
— e 
2tt 



ipi 



G{p). 



(7) 



It can be shown (see for instance [8j) that the posi- 
tivity of the spectral function implies the positivity of 
C(t). In Fig. [3j we plot a numerical Fourier transform of 
the SU(3) gluon propagator for the best parameters de- 
scribed above: this shows a clear violation of positivity. 
We observe that the curve of C (t) is very similar to the 
one of [7|: it is strongly positive for t <1 fm and slightly 
negative beyond. 

Let us now consider the three-dimensional case, where 
the 1-loop calculation can be done explicitly again. The 
details of the calculations and the final expressions will 
be presented in a future publication [24j. We only men- 
tion here that the model ^ is able to account for the 
main features of gluon and ghost propagators found in 
lattice simulations. In Fig. 0] the results of the present 
model with the best fit parameters g — 3.7 \J GeV and 
to = 0.89 GeV for /1 = 1 GeV are compared with d = 3 
simulations performed with the gauge group SU(2) 0]. 
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FIG. 4: Three-dimensional functions for SU(2) gauge group: 
Comparison of present results (plain red curve for jj, = 1 GeV 
and dashed blue curve for fi = 11 GeV) with lattice data of 2] 
(green bars). Top figure: gluon propagator. Bottom figure: 
ghost dressing function. 



We observe that the best fit for gluon and ghost prop- 
agators are not as good as in d — 4. This is probably 
related to the fact that higher loop corrections are not 
very small. It is worth mentioning that the results im- 
prove if one imposes the normalization conditions at a 
larger momentum scale (for n — 11 GeV the best param- 
eters are g = 1.6 \J GeV and m — 0.35 GeV). Such a large 
scheme dependence indicates that higher loop corrections 
give significant contributions. In any case, our calcula- 
tion reproduces the finite IR gluon propagator and ghost 
dressing function. It also reproduces the non-monotonic 
behaviour of the gluon propagator in the IR. An expan- 
sion of the inverse propagator at low momentum leads to 
m 2 - Ng 2 p/64 + 0{p 2 ). 

Within the model, the difference between d = 2 and 
d > 2 that is observed in the lattice (see above) also ap- 
pears natural. In d — 2, we find that the gluon mass 
m and ghost dressing function at zero momentum F(0) 
develop logarithmic IR divergences. Such divergences ex- 
clude the possibility of controlling the one-loop calcula- 
tion as was used above for d > 2. A proper treatment of 
the d = 2 case requires a renormalization group approach 
adapted to the IR regime |24j and goes beyond the scope 
of the present letter. 

The model presented here reproduces surprisingly well 
the d — 4 propagators of pure gluodynamics for SU(2) 
and SU(3) and describes in a simple way the main charac- 



teristics of those propagators in d = 3. The specificities of 
the d = 2 case result from the IR divergences that appear 
in this dimension. Given the technical simplicity of this 
approach, this work opens the door to many subsequent 
applications in Strong Interactions physics (for example, 
the inclusion of quarks, a study of the dependence on the 
gauge fixing, three and four point correlation functions, 
quark- anti-quark static potential). 

Considering the surprisingly good agreement between 
the 1-loop calculation and the simulations, it is tempting 
to think that the model is not just a good phenomeno- 
logical description and one should try to justify the use 
of this action from first principles. The issue of unitarity 
should be also explored. 
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